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Abstract 

We discuss the possibility to probe the QCD critical point during the dynamical black hole formation from a gravitational collapse 
of a massive star, where the temperature and the baryon chemical potential become as high as T ~ 90 MeV and fig ~ 1300 MeV. 
Comparison with the phase boundary in chiral effective models suggests that quark matter is likely to be formed before the horizon 
is formed. Furthermore, the QCD critical point may be probed during the black hole formation. The critical point is found to move 
in the lower temperature direction in asymmetric nuclear matter, and in some of the chiral models it is found to be in the reachable 
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region during the black hole formation processes. 
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1. Introduction 

The critical point (CP) of the Quantum Chromodynamics 
(QCD) |1] may be regarded as a corner stone of the QCD phase 
diagram; the cross over and the first order phase boundaries 
between the hadron and quark-gluon phases are connected by 
CP, then it determines the global structure of the phase diagram. 
Against its importance, the existence, number and location of 
CP are not yet established in theoretical calculations. The sign 
problem makes the lattice Monte Carlo simulation difficult in 
the large baryon-chemical potential (u B ) region J^, the strong 
coupling expansion of lattice QCD Q4£| is not yet reliable to de- 
scribe the real world, and the predictions in effective models 
are spread in the T - ps plane fl. Thus further experimental 
and theoretical developments are necessary to reveal the prop- 
erties of CP. The search for CP in heavy ion collisions is on- 
going in the beam energy scan programs at RHIC and CP 
is one of the most important targets in the forthcoming FAIR 
facility. The most characteristic feature of CP is the divergence 
of the coherence length The phase transition becomes the 
second order, and large fluctuations of the order parameters are 
expected in a volume of the size £ 3 '. On the basis of this idea, 
various signatures of CP have been suggested theoretically U 
However, since the system size and the evolution time are lim- 
ited, it is not an easy task to observe the divergence signature of 
£ in heavy-ion collisions Q]- In addition, if the baryon chem- 
ical potential of CP (pep) is above 500 MeV, CP may not be 
reachable in heavy-ion collisions. Therefore, it is important to 
examine other candidate sites where hot and dense matter is 
formed and CP is reachable. 

A gravitational collapse of a massive star may be one of the 
promising candidates for the CP hunting. It has been argued 
that the transition to quark matter might trigger the second col- 



lapse and bang in supernova explosions |0,[l(i. Recent calcula- 
tion suggests that the QCD phase transition may take place dur- 
ing the core-collapse of a star with M = (10 - 15)Mp when one 
uses the relativistic equation of state (Shen EOS) 11111 combined 
with the bag model EOS at high T or ps for a small bag con- 
stant, B l/4 a 162 MeV The QCD phase transition leads to 
the second shock, and is suggested to give successful supernova 
explosion even in a simulation with spherical symmetry. Non- 
rotating massive stars with mass M > 20M o are expected to 
collapse without supernova exp losions and to form a black hole 
(BH) In Refs. QQQjL the BH formation processes are 
calculated by using the neutrino-radiation hydrodynamical sim- 
ulations in general relativity. In the collapse and bounce stage 
of a 4OM star, the core bounce launches the shock wave, but 
the shock wave stalls due to the collisions with falling matter 
(accretion) and goes down to the surface of the compact object. 
The proto-neutron star is born at center and gradually contracts. 
Because of the accretion, the proto-neutron star mass increases 
rapidly and reaches the critical mass. The dynamical collapse 
occurs again at this point and the BH is formed at 1.3 s af- 
ter the bounce in the case of the Shen EOS 111 111 . If we adopt 
a relativistic EOS with hyperons (Ishizuka EOS) IU6I1 or the 
Lattimer-Swesty EOS 1 1711 . the second collapse becomes more 
rapid (~ 0.7 s after the bounce), while the average neutrino en- 
ergies are lower with Ishizuka EOS 111 51 . The combination of 
the neutrino duration time and the neutrino energy may be used 
as a signal of the hyperon emergence or other new degrees of 
freedom during the BH formation [18]. 

During the BH formation, hot (T ~ 90 MeV) and dense 
(Pb ~ 4po) matter is created. The temperature and density in 
BH formation are significantly higher than those in the model 
explosion calculation of supernova, where the highest temper- 
ature and density are (T,p B ) ~ (21.5 MeV, 0.24 frrT 3 ) fu&. In 
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hotter and denser environment during BH formation compared 
with the supernova explosions, we have a larger possibility of 
creating a new form of matter, such as the dense quark matter. 
Thermodynamical variables at a given time vary as a function 
of radius in a proto-neutron star and form a line in the T - pb 
plane. This thermodynamical line, referred to as the BH forma- 
tion profile in the later discussions, evolves with time and may 
pass through CP and the vicinity. We call here this situation as 
the CP sweep. 

In this Letter, we examine the location of the QCD phase 
boundary and CP in two-flavor chiral effective models at finite 
isospin chemical potential dp = (pd - p u )/2, and discuss the 
possibility of the CP sweep during the BH formation from a 
gravitational collapse of a massive star. First, we compute the 
CP location in the T-ps plane in chiral effective models such as 
the Nambu-Jona Lasinio (NJL) model lfl9ll . the Polyakov loop 
extended Nambu-Jona-Lasino (P-NJL) model E I2H I2H1 , P- 
NJL model with eight-quark interaction 112311 . and the Polyakov 
loop extended quark-meson (POM) model 124 12511 . The bag 
model EOS adopted in Ref. [ 10] is not suited to the present 
purpose, since it does not have CP. In the dynamical BH forma- 
tion, we have abundant neutrinos, and approximate /3 equilib- 
rium including neutrinos are realized inside the neutrino sphere, 
Hn ~Pp = Pe -fiv, while neutrinos are out of equilibrium outside. 
In both cases, it is necessary to take account of finite isospin 
chemical potential 6p as another independent thermodynamical 
variable 12611 . rather than imposing the neutrino-less ji equilib- 
rium condition (dp = p e /2) in order to examine the CP property 
during the BH formation. Recently, P-NJL model with isospin 
chemical potential has been investigated with [27] and with- 
out B28I1 neutrino-less f3 equilibrium conditions. At finite dp, 
we naively expect that CP moves in the lower T direction be- 
cause of the larger rf-quark chemical potential fij = //g/3 + dp. 
Since the matter passes through the high ps and low T region, 
the reduction of the CP temperature Tqp is essential for the CP 
sweep during the BH formation. Next, we compare the results 
of the CP location in the chiral effective models with the evo- 
lution of thermodynamical variables (T,pB,dp) during the BH 
formation obtained in Ref. 1113(1 . It should be noted that we 
compare the results of the CP location in chiral effective mod- 
els and the thermodynamical condition {T,/ub) calculated with 
the hadronic EOS. This comparison is relevant, since the ther- 
mal trajectory should be the same even if we use the combined 
EOS of quark and hadronic matter, as long as the hadronic EOS 
is reproduced at low T and ps in the combined EOS. Finally, 
we discuss the possibility of the CP sweep during the BH for- 
mation from the above comparison. 



2. Polyakov loop extended chiral effective models 



2.1. NJL model 

The Lagrangian density of the two flavor NJL model is given 

by 

Xnjl = q (iyd^ - m ) q + G IT [(qq) 2 + (qijsrq) 2 ] 
- G p [(qy^rqf + (qiysyrq) 2 ] - G u [(qy^q) 2 + iqiy^q) 2 ] , 

(1) 

where q denotes a quark field with Dirac, color and flavor in- 
dices; t is the Pauli matrices in the flavor space. In what fol- 
lows, we take G p = G M = Gy, which amounts to take <d and p 
mesons degenerate in the vacuum. We can fix the value of Gs 
by fitting the well known properties of the QCD vacuum. We 
will take the ratio Gy/Gs as a free parameter. 

We are interested in unbalanced populations of u and d 
quarks, not in the neutrino-less /^-equilibrium. For this reason, 
we introduce two independent chemical potentials for u and d 
quarks. At the mean field level, the effect of the vector inter- 
action is to shift the quark chemical potentials: the flavor sin- 
glet interaction gives a contribution proportional to p u + pd\ on 
the other hand, the flavor triplet interaction gives a contribution 
proportional to the isospin density, p u - p c i- Keeping this into 
account, and for later convenience, we define 



p n = p-dp- AG v p u , p d = p + 5p- 4G v pd 



(2) 



where p and dp represent chemical potentials conjugated to the 
total quark number density and to isospin density, respectively. 

The one-loop thermodynamic potential can be represented as 
the sum of the the vacuum Qo and the thermal (finite tempera- 
ture and finite chemical) Q7 contributions, 



Qnjl — Qo + Qr. 

r2 



Or = -2riV c X/^log(l + ^)(l +e -^) 
-2G v (p 2 l +p 2 l ), 



(3) 
(4) 



&{ =E p ±p f , 



(5) 
(6) 



where E p = ^p 2 + M 2 with M = m + 2, and 2 = -2G a (iqq) 
corresponds to the mean field quark self-energy. In Eq. (|4), 
we have treated the ultraviolet divergence of the vacuum en- 
ergy by the use of a smooth regulating function, F(p 2 , A) = 
[1 + (p 2 ) 5 /A 10 ] _1 . This is done just for numerical convenience; 
quantitatively, we have checked that the results are consis- 
tent with the more common hard cutoff regularization scheme, 
within a few percent. The thermal part Q7- is a finite contribu- 
tion which does not need any regularization. 



In this section, we summarize the chiral effective models, the 
NJL, P-NJL, and PQM models, which we use in computing the 
CP location in the T - ps plane. 



2.2. The P-NJL model 

The P-NJL model Lagrangian density is still specified by 
Eq. ([TJ, with the derivative replaced by a covariant one: d M — > 
Dft = d M - iA M . Here, A p is a temporal, static and homogeneous 
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background gluon field related to the Polyakov loop P, whose 
expectation value is computed self-consistently. The one-loop 
thermodynamic potential in Eq. © is replaced with 112 ill 



Q P _ NJL =£lo + £l T + t((P,P, T), 
a(T): 



U[P,P,T] =r 4 |-^-PP + b(T)a\H(P,P)^ , 
Pexpli J ^tA 4 | 



1 



P =— Tr 



= -Tie 

3 



(7) 
(8) 

(9) 



where U(P, P, T) is the Polyakov loop effective potential [22] 
with H(P, P) = 1 - 6PP + 4(P 3 + P 3 ) - 3{PP) 2 , a(T) = 
aQ + a l (T /T) + a 2 (To/T) 2 ,andb(T) = biiTo/T) 3 . We adopt the 
Polyakov gauge, where <p is specified by = 03/I3 + <p%A%. The 
standard choice of the parameters reads II22|| ao — 3.51 , a\ — 
-2A1 , «2 = 15.2, £>3 = -1.75. The parameter To in 
Eq. ([8]l sets the deconfinement scale in the pure gauge theory, 
i.e. T c = 270 MeV. 

The thermal part of the thermodynamic potential is now 
given by 12 ill 



Cl T = - 2G V (pi + P$-2tZJ^ log (f{f£) , (10) 



F{ =1 + 3Pe~ fss+ + 3Pe- 2l3E * + e' 3 ^ , 
Ft =1 + 3Pe-^'- + 3Pe- 2 > 3E - + e^- . 



(11) 
(12) 



In Eq. ( fTTT ). the addenda on the r.h.s. correspond to the ther- 
mal contribution of zero, one, two and three quark states, re- 
spectively. Analogously, Eq. (fT2l is the thermal contribution of 
antiquarks. 

2.3. P-NJL model with eight-quark interaction (P-NJL%) 



With eight-quark interaction, the quark Lagrangian is Q23 1 

,2 



-C =£njl + Go-8 [(qqf + (qijsTq) 2 ] 



(13) 



We do not include here the eight-quark interaction in the vector 
channel. The one loop thermodynamic potential is still given 
by Eq. (0, with 



= 3G^ 4 V_ _ y r _fP F{ 2 m (H) 

16G4 4G (r ^ J (2tt) 3 ^ ' '" 



and E p = y/p T TM 2 , with M = m + 2 + I 3 (G fr8 /2G^.). For 
simplicity, we call this model P-NJLg. 

P-NJLg parameters (A, G^, Go-%, mo) are determined to fix f n , 
m„, m lT , and (uu), and as a result, the chiral and deconfine- 
ment transition temperatures are found to roughly agree 023 1. 
The eight-quark interaction is found to make chiral cross over 
sharper at fi — 0; fi makes this cross over stronger, and eventu- 
ally leads to the first order transition. Hence eight-quark inter- 
action causes chemical potential of CP to be smaller. 



2.4. PQM model 

The Lagrangian density of the PQM model is given by 

£ =q [i^Dp - g(cr + iy 5 T ■ n) - gvfiUjx + t ■ Rp)] q 

+ \{dp(r) 2 + l -{dpnf - Via, n) - U(P, P, T) 

Aojp V of v - l -Rp V ■ R" v + l -m 2 ,{u>puf + Rp ■ /?") . (15) 

The mesonic potential is U (o~, n) — A (cr 2 + n 2 - v 2 ) 14 -ha, 
and LOp Y and Rp V are the field tensors of the a> and p mesons. 
We use the same Polyakov loop effective potential as that in 
the P-NJL model, Eq. ([8]). In the mean field approximation, the 
thermodynamic potential is given by 



n PQM =14(P, P, T) + U(cr, n = 0) + Q + O r 



(16) 



/d 3 
—^Ep9(A 2 -p 2 ), (17) 
\2ny 

where Q.q corresponds to the regularized fermion vacuum en- 
ergy with M = gcr, and D.j is still given by Eq. ( fTOb with 
Gy — g 2 J2m 2 . While the PQM model is renormalizable and 
an elegant procedure of dimensional renormalization is feasi- 
ble 125J, it is enough to cut large momenta by a hard cutoff for 
our purposes. 



2.5. Model parameterization 

In this study, we fix fi and 6fi, and compute p„, pd, the chi- 
ral condensate and the Polyakov loop expectation values self- 
consistently, requiring the stationary condition of the thermo- 
dynamic potential. In the case of the NJL model and the P- 
NJL model without eight-quark interaction, the parameters Go-, 
A and mo are chosen in order to reproduce the QCD vacuum 
properties (uu) = (-250 MeV) 3 , /„ = 92.4 MeV and m n = 139 
MeV. They are given as A = 618.98 MeV, G a = 2. 05 /A 2 , 
and mo = 5.28 MeV. The parameter To in the Polyakov loop 
effective potential is taken to be To = 210 MeV. With this pa- 
rameter choice, the constituent quark mass in the vacuum is 
M ^ 340 MeV. In P-NJLg, we use the parameterization in [23], 
A = 631.5 MeV, G a = 1.864/A 2 , G lT& = 11. 435/ A 8 , and 
mo — 5.5 MeV, which give the vacuum constituent quark mass 
M a; 353 MeV. For Gy, estimates exist based on perturbative 
one gluon exchange i29ll . r = GvlG a — 0.5; on instanton-anti- 
instanton molecule model i30ll . r — 0.25; and an interpolation 
is obtained by a lit of the Lattice data with the P-NJLs model 
in Bill , r = 1. We here treat r as a free parameter, and compare 
the results with r — and r = 0.2. 

The parameters of the PQM model, v, A, g and h, are fixed 
to reproduce some vacuum properties: the chiral condensate in 
the vacuum, cr — f„ — 92 A MeV; the vacuum pion mass, m 2 = 
h/f K = (139MeV) 2 ; the constituent vacuum quark mass, M = 
gf n = 335MeV; the sigma mass, given by m 2 = d 2 Q/<9<x 2 = 
(700 MeV) 2 . In this article, we use the following parameter 
set: A = 600MeV, v 2 = -(617.68MeV) 2 , and A = 2.7255. 
We use the same Polyakov loop effective potential as that in the 
P-NJL model. The vector meson mass is chosen to be m v = 
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770 MeV. The vector coupling is treated as a free parameter, 
and we compare the results with r = g v /g = and 0.2 in the 
later discussions. 

P-NJL models with SU(3) / are known to support hybrid neu- 
tron star mass of 1 .74-2. 12 M Q when combined with a hadronic 
EOS at low densities lL32ll . By comparison, an SU(3)/ version 
of PQM model connected with a hadronic EOS can support the 
hybrid neutron star mass of 2.0M o A3 311 . Generally speaking, 
the EOS becomes softer with larger degrees of freedom, then 
SU(2)/ versions of these models are expected to give stiffer 
EOS than in SU(3)/. Therefore SU(2)y chiral effective models 
adopted in this work would result in larger maximum masses of 
neutron stars, which are consistent with the recently observed 
1.97 + O.O4M neutron star Q. 



Table 1: Location of CP, the transition chemical potential at T = (fi c ), and 
the type of the transition to quark matter during the BH formation. All T and /i 
values are given in the unit of MeV. 



3. Critical point location and its sweeping 
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Figure 1: Phase diagram in (T,fiB,S^t) space. First order phase boundaries 
{T,/ib) calculated with PQM are shown for several values of the isospin chem- 
ical potential, S/i. We also show the BH formation profile, (thermodynamical 
profile (T,/jb,S/i) during the BH formation) at t = 0.5, 1.0, 1.344 sec after the 
bounce (double lines). 



3.1. Critical point and phase boundary in asymmetric matter 

In Fig. Q] we show the isospin chemical potential depen- 
dence of the first order phase boundary and the critical point 
in the PQM model. We find a trend that the first order phase 
boundary shrinks at finite isospin chemical potential. Transi- 
tion temperature at a given baryon chemical potential ps = 3p 
decreases, and the transition chemical potential fx c at T — 
also decreases. We do not consider here the pion condensed 
phase, because the s-wave pion condensation will not be real- 
ized when we include the s-wave ttN repulsion 13511 , The CP 
location is sensitive to dp. Compared with the results in sym- 
metric matter, Tqp becomes smaller at finite dp and reaches zero 
at dp = dp c =s (50 - 80) MeV. The downward shift of Tqp may 
be understood from the density shift. At low T and without 
the vector interaction, the quark density is proportional to p 3 , 
Pu,d ^ Qi + dpf . Then the sum of u and d quark density in- 
creases when dp is finite, and it simulates higher fi, where the 
transition temperature is lower. 
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In TableQ] we summarize the CP location (T,pb) for several 
values of dp and r = G v /G lT in NJL, P-NJL, P-NJL 8 , and PQM 
models. In P-NJL 8 , our results at r — are in agreement with 
those in Ref. lEill . The transition chemical potential at T — 
is in the range of 1000 MeV < p c < 1110 MeV. p c is sensi- 
tive to the details of the interaction, especially to the strength 
of the vector interaction. The temporal component of the vec- 
tor potential shifts the chemical potential effectively as already 
introduced in Eq. ©. In the momentum integral, we find the 
effective chemical potential pj = p + 5(1 — Vf appears, where 
Vf = 4-GvPf represents the vector potential for quarks. The re- 
pulsive vector potential reduces the effects of the chemical po- 
tential and consequently leads to an upward shift of p c by about 
10-15 MeV at r — 0.2. When we increase the vector coupling 
from r — to r — 0.2 in NJL, the first order transition boundary 
is shifted upward in p and 7cp is reduced from 50 MeV to 22 
MeV. At larger vector coupling, the first order phase boundary 
disappears, and the QCD phase transition becomes the cross 
over at any p. This trend also applies to the P-NJL and PQM 
models; the phase boundary is shifted in the larger p direction 
and shrinks in the T direction with finite vector interaction. In 
P-NJLg, the first order phase transition is robust and survives 
with larger vector interaction such as r — 0.8, while the effects 
of the vector interaction is qualitatively the same. 

3.2. BH formation profile 



In Fig. [2] we show the BH formation profile {T,ps, 5p) 1 1311 
calculated by using the Shen EOS at t = 0.5, 1.0 and 1.344 sec 
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Figure 2: The BH formation profile, (T.fig. S/i), as a function of the radius. We also show the baryon density (pg) and the electron fraction ( Y e = Pe/ps)- Results are 
shown for the gravitational collapse of a 40 M Q star at t = 0.5 sec (dotted lines), 1.0 sec (dashed lines), and 1.344 sec (solid lines, just before the horizon formation). 



after the bounce during the BH formation from a 40 M star in 
the proto-neutron star core, where the mass coordinate from the 
center is M < 1.6M Q . The time t = 1.344 sec is just before 
the horizon formation. From the outer to the inner region of 
the proto-neutron star, T first increases from T ~ 10 MeV to 
T ~ (50 - 90) MeV in the middle heated region, and decreases 
again inside. The central density grows from ps ~ po at bounce 
to 2po, 2.5po and 4po at 0.5, 1.0 and 1.344 sec, respectively. 
The charge to baryon ratio (Y e ) is less than 0.3 inside the proto- 
neutron star IU4II . The isospin chemical potential is found to 
be 50 - 130 MeV in the inner region. The baryon chemical 
potential ps is found to go over 1300 MeV in the central region 
just before the horizon formation at t — 1 .344 sec. 

The maximum chemical potential is much larger than the 
A(l 1 15) mass, and hyperons are expected to emerge. Actually, 
hyperons are formed abundantly when we use Ishizuka EOS 
including hyperons ifTsll . while the proto-neutron star collapses 
earlier and the maximum ps (~ 1 100 MeV) is lower. 

It should be noted that the maximum mass of neutron stars 
is 2.17 M in Shen EOS, which does not contradict to the ob- 
servation of 1.97M neutron star 113411 . We need care in dis- 
cussing the results of Ishizuka EOS, since the maximum mass 
is 1.63M and it cannot support the 1.97M mass neutron star. 
It would be important to discuss the hyperon admixture with 
more repulsive hyperon potentials at high densities, which sup- 
port heavier neutron stars. 

3.3. Possibility of critical point sweep 

We shall now compare the CP location and the phase bound- 
ary with the BH formation profile. In Fig. [3] we compare 
the phase boundaries and the CP location in NJL (left-top), P- 
NJL (right-top), P-NJL 8 (left-bottom) and PQM (right-bottom) 
models with the BH formation profile, (T,ps). As already 
mentioned, the transition chemical potential is in the range of 
1000 MeV < fx c < 1110 MeV in symmetric nuclear matter, and 
it decreases at finite 6p. During the BH formation, the baryon 
chemical potential reaches around 1000, 1100 and 1300 MeV 
in the central region of the proto-neutron star at t — 0.5, 1 .0 and 
1.344 sec, respectively. This comparison suggests that quark 
matter would be formed between t - 0.5 and 1.0 sec in the 
central region of the proto-neutron star in most of the models 
considered here. 

Since the CP location has strong dependence on model and 
parameter, there are three possible types in the transition to the 



quark matter during the evolution of matter toward the BH for- 
mation; the first order transition, the cross over transition, and 
the CP sweep. In P-NJL^ (with and without vector interaction), 
Tcp is relatively high even in asymmetric matter, then the mat- 
ter experiences the first order phase transition, and CP is not 
reached. In NJL and P-NJL with vector interaction, Tqp de- 
creases in asymmetric matter, and CP already disappears in the 
central region at t = 1.0 sec, where the isospin chemical po- 
tential is large, dp ~ 70 MeV. In this case, the BH forma- 
tion profiile evolves above CP, and the cross over transition to 
quark matter will proceed without going through the first or- 
der boundary. In NJL and P-NJL models without vector inter- 
action and PQM with and without vector interaction, the BH 
formation profile goes through CP from below, as shown in the 
double line in Fig. \T\ CP in symmetric matter is above the BH 
formation profile at t — 1.0 sec, while CP in asymmetric mat- 
ter is below the line at t — 1.344 sec. Since the matter in the 
central region is highly asymmetric {dp = (50 - 130) MeV) at 
t = 1.344 sec, some part of the off-center BH forming matter 
would go through CP between t — 1 and 1.344 sec, i.e. CP is 
swept. 

CP sweep will take place in the off-center inner core. As can 
be seen in Fig. [2] the highest density and baryon chemical po- 
tential is realized at center, while temperature becomes highest 
off-center due to the shock wave. In the PQM (r - 0) case, 
as a typical example of CP sweep, the quark matter formation 
takes place first at center with T < Tqp, the quark matter re- 
gion grows, and the CP sweep takes place when the boundary 
temperature reaches 7cp. We find from Fig. [3] that CP will be 
swept at a little inner region than the highest T radius, and this 
corresponds to R — 6 - 8 km as seen in Fig. [2] 

It is also interesting to discuss whether the condition in 
core-collapse supernovae may cross the phase boundaries to 
the quark matter phase. We examined the numerical result of 
the gravitational collapse and bounce from the massive star of 
15M with Shen-EOS l36ll . The density and temperature at the 
core bounce are (p B , T) = (0.2fnT 3 , 15 MeV) at center. The 
baryon chemical potential and temperature do not reach the 
phase transition region in this case. They do not rise enough 
during the stage of the stalled shock wave for 200 ms after the 
core bounce. Therefore, it is not likely to have the CP sweep in 
dense matter around the core bounce for successful explosions. 
The situation is similar in the numerical result of adiabatic col- 
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Figure 3: Critical point, phase boundary and the BH formation profile. Critical point in symmetric (circles) and asymmetric (squares) matter and the first order 
phase transition boundaries (solid lines) in chiral effective models are compared with the BH formation profile (T,fig) at t = 0.5, 1.0 and 1.344 sec (dashed lines). 
Top-left, top-right, bottom-left and bottom-right panels show the results in the NJL, P-NJL, P-NJLs, and PQM models, respectively. We show the results without 
(r = 0) and with (r > 0) vector interaction. In P-NJL with vector interaction and PQM, we also show the CP trajectory (dotted lines). Arrows around the BH 
formation profile show time evolution, and arrows around the critical points show the isospin chemical potential dependence. 



lapse and bounce of the 15M star using Ishizuka' hyperonic 
EOS as discussed in 111 611 . However, there are still possibili- 
ties to have the quark matter and the CP sweep in the late stage 
of core-collapse supernovae. In the thermal evolution of the 
proto-neutron star for ~ 20 sec after the birth in supernovae, 
the density and temperature increase further due to the contrac- 
tion during the thermal evolution. In fact, the condition may 
become closer to that of CP in the long term evolution after the 
bounce [36]. In the massive proto-neutron star cases ll3~7il . hy- 
perons and quarks appear in the late stage due to the extreme 
conditions, where they reach the phase boundaries. Hence, the 
CP sweep may be realized in the evolution toward the black 
hole formation in the very massive proto-neutron stars, while 
the explosion goes on. It is also possible to have high density 
and temperature in the transition region to quark matter at the 
second bounce of the core collapse supernova when one adopts 
an EOS having a low transition baryon chemical potential 



lopts 



There are two comments in order, (a) While we have 
compared the CP location and the BH formation profile in 
(T,p B , Sp), the relation between the densities and chemical po- 
tentials is highly model dependent and the baryon densities in 
the quark and hadronic models considered here do not gen- 
erally coincide; For example, in NJL and P-NJL, the baryon 
densities on CP are found to be p B = 0.30 and 0.29 fnT 3 at 
dp = 65 MeV, respectively, while p B = 0.36 and 0.31 fm~ 3 
in Shen EOS at the corresponding (T,p B ,6p). In PQM with 
r = and 0.2, p B = 0.25 and 0.23 fnT 3 at dp = 60 MeV, 
respectively, while p B = 0.33 and 0.34 fm 3 in Shen EOS on 
these points. Thus the baryon densities in chiral effective and 
hadronic models have 10-50% differences in density at CP, and 
the results shown in this article would be modified when we per- 
form a consistent calculation using a hybrid supernova matter 



EOS which smoothly connect quark and hadronic matter EOSs. 
However, the trajectories span a wide region in (T,p B ) which 
are in a range expected for the phase transition lines and CP, 
then we expect the CP sweep or at least the transition to quark 
matter even in a consistent calculation of the black hole forma- 
tion, (b) The symmetry energy in RMF models has a strong 
dependence in p B and may overestimate the effects at high den- 
sities. When one adopts a weak p B dependence of the symmetry 
energy at high densities, smaller isospin chemical potential re- 
sults. In this case, the cross over transition during BH formation 
becomes less probable, since Tqp is kept to be high at small bp. 

4. Summary 

In this Letter, we have discussed the possibility of the QCD 
phase transition to quark matter and the critical point (CP) 
sweep during the dynamical black hole (BH) formation. We 
have compared the phase boundary and CP in chiral effective 
models with the BH formation profile, thermodynamical vari- 
ables (T,p B ) calculated in the neutrino-radiation hydrodynam- 
ics. For this comparison, it is necessary to consider asymmetric 
matter at finite isospin chemical potential, dp = {p„-p p )/2 + 0. 
The isospin chemical potential is found to reduce the tempera- 
ture of the critical point Tcp, then we have a larger possibility 
of the CP sweep or the cross over transition to quark matter. 

In the chiral effective models considered here, with and with- 
out the vector interaction, the transition chemical potential at 
T = is found to be in the range of p c = (1000 - 1 1 10) MeV in 
symmetric matter, and p c decreases at finite dp. We can com- 
pare these values with the highest baryon chemical potential 
realized during the BH formation, p B = 1300 MeV. Then if the 
baryon chemical potential is larger than the QCD phase transi- 
tion chemical potential, quark matter will be formed during the 
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BH formation. In order to conclude, however, it is necessary to 
examine with the EOS which includes both baryonic and quark 
degrees of freedom. The CP location is sensitive to the mod- 
els and parameters. We have found that there are three types of 
possibilities of the transition to quark matter in the BH forma- 
tion process. When the thermodynamical trajectory go below, 
above, or through CP in asymmetric matter, the QCD phase 
transition proceeds via the first order transition, the cros s over 
transition, or the CP sweep. When CP is swept, the density fluc- 
tuation should grow, and we can expect various effects around 
CP such as the baryon density fluctuations, focusing of ther- 
modynamical trajectories, disappearance of the sound mode, as 
discussed for the beam energy scan at RHIC 10]. 

It is a big challenge to construct an EOS which is applica- 
ble to the dynamical simulation of the core-collapse processes, 
contains both the hadronic and quark degrees freedom, and in- 
cludes CP of QCD. There is an attempt to include both quark 
and baryon contributions based on the PQM ll38ll . In that work, 
the vacuum quark contribution is ignored and the QCD phase 
transition becomes the first order in the two-flavor chiral limit. 
It may be also necessary to consider the effects of inhomoge- 
neous phases 113911 . which may emerge around the first order 
phase boundary. 
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